Extended quantum f/(l) -liquid phase in a three-dimensional quantum dimer model 
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Recently, quantum dimer models have attracted a great deal of interest as a paradigm for the study of exotic 
quantum phases. Much of this excitement has centred on the claim that a certain class of quantum dimer model 
can support a quantum f/(l) liquid phase with deconfined fractional excitations in three dimensions. These 
fractional monomer excitations are quantum analogues of the magnetic monopoles found in spin ice. In this 
article we use extensive quantum Monte Carlo simulations to establish the ground-state phase diagram of the 
quantum dimer model on the three-dimensional diamond lattice as a function of the ratio p, of the potential 
to kinetic energy terms in the Hamiltonian. We find that, for jic = 0.75 ± 0.02, the model undergoes a first- 
order quantum phase transition from an ordered "R-state" into an extended quantum (7(1) liquid phase, which 
terminates in a quantum critical "RK point" for ^ = 1. This confirms the published field-theoretical scenario. 
We present detailed evidence for the existence of the (7(1) liquid phase, and indirect evidence for the existence 
of its photon and monopole excitations. Simulations are benchmarked against a variety of exact and perturbative 
results, and a comparison is made of different variational wave functions. We also explore the ergodicity of the 
quantum dimer model on a diamond lattice within a given flux sector, identifying a new conserved quantity 
related to transition graphs of dimer configurations. These results complete and extend the analysis previously 
published in [O. Sikora et al. Phys. Rev. Lett. 103, 247001 (2009)]. 

PACS numbers: 75.1().Jm 75.10.Kt, ll.15.Ha, 71.10.Hf 



I. INTRODUCTION 

Dimer models, which describe the myriad possible con- 
figurations of hard-core objects on bonds, have long been a 
touch-stone of statistical mechanicsi*^. Quantum dimer mod- 
els (QDM's), in which dimers are allowed to resonate between 
different degenerate configurations, were first introduced to 
describe antiferromagnetic correlations in high- Tc supercon- 
ductorsi. It has since been realized that QDM's aiise naturally 
as effective models of many different condensed matter sys- 
tems, and provide a concrete realizations of several classes 
of lattice gauge theories. As such, they have become central 
to the theoretical search for new quantum phases and excita- 
tions. A key question in this context is when, if ever, a QDM 
can support a liquid ground state? 

The answer to this question depends on lattice dimension 
and topology. In two dimensions (2D), the situation is rela- 
tively well-understoodi, and a liquid ground state is known 
to exist in the QDM on the non-bipartite triangular-5. and 
kagome^ lattices. This liquid is gapped and has deconfined 
fractional excitations, which are vortices of an underlying Z2 
(Ising) gauge theory. Meanwhile, for 2D bipartite lattices, 
the underlying gauge theory has a C/(l) character, akin to 
electromagnetism^'*^. In this case a liquid state is found only 
at a single, critical, "Rokhsar-Kivelson" (RK) point^. Away 
from this, the system crystallizes into phases with broken lat- 
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FIG. 1 : (Color online) Phase diagram of the quantum dimer model 
on a diamond lattice, as a function of the ratio p of potential to kinetic 
energy, following**''^''"*. 



tice symmetries. For the QDM on a square lattice, these are 
columnar dimer and resonating plaquette phasesMi. A very 
similar phase diagram is found for the QDM on a honeycomb 
latticaii. 

Much less is known about QDM's in three dimensions (3D). 
However here field theoretical arguments suggest that a new, 
extended C/(1)-Iiquid phase with gapless photon-like excita- 
tions might "grow" out of the RK point for a 3D QDM on 
a bipartite lattics^ii^ (see Fig. [T]). Similar claims have been 
made for the closely related quantum loop model in three di- 
mensions'-. 

In a recent Letter—, we presented numerical evidence for 
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the existence of an extended quantum U{1) liquid phase 
with deconfined fractional excitations in a three-dimensional 
QDM, taken directly from its microscopic Hamiltonian. We 
considered the specific example of a QDM on a diamond 
lattice, which may serve as an effective model for the frus- 
tration found in various spinel compound a'^i'^ . From zero- 
temperature GFMC simulations of clusters of up to 2000 sites 
we were able to demonstrate the existence of a quantum U{1)- 
liquid phase for a range of parameters 0.77 ± 0.02 < /i < 1 
bordering on the RK point [cf. Fig. [II- These results comple- 
ment a recent finite-temperature quantum Monte Carlo study 
of a microscopic model of interacting bosons which can be 
mapped onto an effective quantum loop modeU^. 

In this article we document the evidence for these claims, 
taken from explicit calculations of (i) the order parameter of 
the competing (dimer)-ordered phase [cf. Fig. |2l, (ii) the 
"string tension" associated with fractional monomer excita- 
tions and (iii) the finite-size scaling of a characteristic set of 
energy gaps associated with the U(l) liquid phase. We also 
explore the internal consistency of these calculations, cross- 
checking simulations against exact diagonalization of small 
clusters and perturbative expansions about exactly soluble 
limits of the model. 

We pay particular attention to the thorny question of the er- 
godicity of GFMC simulations in the presence of hidden quan- 
tum numbers. In the 2D QDM's studied previously, the quan- 
tum numbers are well understood, and it is clear which dimer 
configurations are connected by the QDM Hamiltonian. This 
is not, however, the case for the 3D QDM considered in this 
article. Despite the fact that we can define a "magnetic" flux 
(j) which is a good quantum number — in a manner similar to 
the square-lattice QDM in 2D — the QDM Hamiltonian on 
a diamond lattice connects only sub-sets of dimer configura- 
tions within of each flux sector. This means that the QDM on 
diamond lattice is not ergodic within a given flux sector, and 
great care must be taken to ensure that simulations accurately 
represent the true ground state for each value of flux 0. We 
shed some light on this issue by identifying a new, hidden, 
quantum number conserved by the QDM on diamond lattice. 

The paper is structured as follows : 

In Section HI] we introduce the model and summarize the 
properties of the proposed ?7(l)-liquid phase. We demon- 
strate a way of constructing dimer configurations with differ- 
ent values of magnetic flux, and define an order parameter for 
the ordered "R-state" found for negative values of ^. In Sec- 
tion [III] we use Green's function Monte Carlo simulations to 
obtain explicit predictions for the fi dependence of the order 
parameter, monomer string tension and ground state energy 
as a function of magnetic flux for a range of different cluster 
sizes and geometries. We confirm the conjectured form of the 
phase diagram Fig.[T] and use careful finite size scaUng to ob- 
tain an improved estimate of the location of the phase transi- 
tion from ordered to liquid phases in the thermodynamic limit 
as He = 0.75 ± 0.02. In Section HVl we explore the ergodic- 
ity of our simulations within a given flux sector and confirm 
that the techniques used give reliable estimates of ground state 
properties. We also identify an additional Z2 quantum number 
which divides each flux sector into two distinct parts. In Sec- 



tion|V]we summarise our results and briefly discuss remaining 
open issues. 

The paper concludes with some technical appendices. In 
Appendix |A] we discuss the generalization to a fermionic 
QDM and show that it is equivalent to the bosonic case in low- 
est order. In Appendix IB] we exhibit a Landau theory for the 
ordered, R-state. In Appendix |C] we compare different vari- 
ational wave functions, and use these to investigate the tran- 
sition from ordered to liquid ground states. In Appendix |D] 
we benchmark our Green's function Monte Carlo simulations 
against perturbation theory about exactly soluble limits of the 
model. 



II. THE QUANTUM DIMER MODEL AND WHAT YOU 
NEED TO KNOW ABOUT IT 

A. The model and its derivation 

In the spirit of Ref . S we study the quantum dimer model 
(QDM) 

Hqdm = -.g/Cf + fij\ff (1) 
where the matrix elements of 

-i:(|0)(0|+|o>(0|) p) 

connect different, degenerate dimer configurations by cycli- 
cally permuting dimers on six-bond "flippable" plaquettes of 
alternating filled and empty links within the diamond lattice 
[cf. Fig. H, while 

= i:(|0)(0|+|0)(0|) <3) 

counts the number of such flippable plaquettes. In what fol- 
lows we will consider g > 0, with — 00 < /i < g. 

This QDM Hamiltonian acts on the « 1.3^/^ different 
dimer configurations which cover the links the A^-bond di- 
amond lattice'^. This Hilbert space breaks up into differ- 
ent subsectors Ac, defined by the sets of dimer configurations 
which are connected by the matrix elements of /Cf . The num- 
ber and nature of these subsectors remains an open problem, 
which will be discussed at length in Section lTVl of this article. 

Crucially, however, all off-diagonal matrix elements of 
Eq. ([T]) are zero or negative. By the Frobenius-Perron the- 
orem, the lowest energy state of the QDM within any given 
subsector of the Hilbert space Ac must then be a nodeless su- 
perposition of dimer configurations 

I V^o)a. = ^ac|c) (4) 
ceAc 

with real, positive coefficients Uc- As a result. Quantum Monte 
Carlo simulations of Eq. ([T]i are not subject to a sign problem. 
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Exactly at the RK point = g, the QDM Hamiltonian can 
be written as a sum of projectors 







0)-|0) (0|-(0 



(5) 



In this case the ground state in any subsector of the Hilbert 
space Ac is the zero-energy eigenstate given by the equally 
weighted sum over all dimer configurations within that sub- 
sector 



(6) 



cSAc 




(b) 



The diamond lattice is the medial lattice of the pyrochlore 
lattice, formed by connecting the centres of all tetrahedra (see 
Fig. |2]i, and and these constrained states are in exact one-to- 
one correspondence with the dimer coverings of a diamond 
lattice. The lowest-order process in degenerate perturbation 
theory which connects these states is a cyclic exchange of 
bosons around the smallest, hexagonal plaquette of the py- 
rochlore lattice. This occurs at third order in the hopping of 
bosons, giving g = 12t^/V^. 

Another physically motivated model which reduces to the 
quantum dimer model on the diamond lattice is the easy-axis 
quantum antiferromagnet on a pyrochlore lattice 



(8) 



at half-magnetization (5^) = S/2. For 5=1/2 this is ex- 
actly equivalent to the hard-core boson problem considered 
above, with each occupied site corresponding to a "up" spin, 
and each empty site to a "down" spin, with the caveat that 
the kinetic energy term J^y now has the opposite sign, and so 
g = -12JIJJI 

At first sight this might seem to imply that the XX Z model 
Eq. ([Hi with .g < has a different ground state from the a.t-V 
model with g > Eq. (|7]). However the global sign of g can 
be changed by the simple unitary transformation 



|c) exp[i7r7VA/2]|c) 



(9) 



where | c) is an arbitrary dimer configuration and Aa counts 
the number of dimers which occupy the subset of diamond 
lattice bonds A shown in Fig. [3] The cyclic permutation of 
dimers on any six-bond flippable plaquette changes A/a by 
two, and so the unitary transformation Eq. (|9]l maps g — > —g. 
Moreover, exactly the same quantum dimer model Eq. ([l]) can 
be derived from a fermionic t-V model, i.e. Eq. ^ with Fermi 
operators substituted for the bosonic ones. This mapping is 
described in Appendix lAl 



FIG. 2: (Color online) (a) Cubic unit cell of the diamond lattice, with 
bonds occupied by dimers shown by thick red links. The dimer con- 
figuration shown is the ordered "R-state" found for negative /i. This 
contains a "flippable" hexagonal plaquette, which is shaded blue, (b) 
Equivalent cubic unit cell of the pyrochlore lattice, with sites occu- 
pied by hard-core bosons shown in red. The boson configuration 
shown is the same ordered "R-state". 

The QDM Eq. ([T]i can be derived from various microscopic 
models. One example is a hard-core boson model on a py- 
rochlore lattice formed of corner sharing tetrahedra 

T^tv = -t ^{h\bj + H.c.) + VY^ n,n.j (7) 

at quarter filling, in the limit of large nearest-neighbor inter- 
actions V ~:^t.\n this case, the low energy configurations are 
those without bosons on neighbouring sites. This imposes the 
constraint of placing exactly one boson in each tetrahedron. 




FIG. 3: (Color online) The subset of diamond lattice bonds A (thick, 
blue, shaded bonds) used in the unitary transformation Eq. (|9j which 
maps g — >■ —g in Eq. (T}. 

An effective Hamiltonian of the QDM form Eq. (HJ can also 
be derived from Eq. ^ for larger values of spin. For S" = 3/2, 
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where the model might be of relevance to Cr spinels, this was 
accomplished by Bergman et al.—, who found ji/g w —7. A 
two-dimensional analogue of Eq. ([B might also be realized 
using cold atoms'"^. Finally, quantum dimer models also arise 
as effective models of spin-orbital models^. 

Since the goal of this paper is to determine the ground state 
phases of Eq. ([l]), and not to relate it to any particular physical 
system, in what follows we set g = 1, and study the model as 
a function of the single adjustable parameter ii. For — )• — cx) 
the potential energy dominates, and the ground state of Eq. ([1) 
is the set of dimer configurations which maximise the number 
of "flippable" hexagons, with one out of four hexagons being 
flippable. This is the so-called R-statei^, illustrated in Fig.|2l 
which has a cubic unit cell and is eight-fold degenerate. 

At the RK point /i = 1, the ground state is the equally 
weighted sum of all possible dimer configurations?, while for 
/i > 1, ground states are the many "isolated states" which 
contain no flippable hexagons, and so have no kinetic energy. 
The interest of the QDM therefore rests in the parameter range 
— 1 < /^j < 1, for which kinetic and potential energy compete 
on an equal footing. It is this range of parameters which we 
consider below. 



B. Effective electrodynamics 

For purpose of comparison with numerics, we now briefly 
review the field theory arguments underlying the proposed 
U{1) liquid phase. The defining property of a dimer model 
is that every lattice site is touched by exactly one dimer The 
diamond lattice is bipartite — we can divide it in two sublat- 
tices, A and B, such that every every link (bond) of the lattice 
connects an A sublattice site with a B sublattice site. We can 
use this property to associate a magnetic flux vector B with 
each bond of the lattice. Bonds not occupied by a dimer are 
assigned one (arbitrary) unit of magnetic flux, directed from 
B to A. Where a dimer is present, we assign three units of 
flux to that bond, directed from A to B^. 

In this language, the constraint that every lattice site is 
touched by one dimer becomes the condition 



V-B = 0. 
We resolve this constraint by writing 
B = V X A, 



(10) 



(11) 



where the gauge field A is defined on the sites of the orig- 
inal diamond lattice while the fictitious magnetic field B is 
defined on its bonds. We chose to work in the Coulomb gauge 
V ■ A = 0. In the absence of quantum effects, the statistical 
description of the dimer model reduces to a relatively straight- 
forward problem in magnetostatics, and the ground state of the 
dimer model is found to be classical [/(I) -liquid. Correlations 
between dimers have a dipolar form in real space, and the con- 
straint V • B = manifests itself as a set of "pinch points" in 
the structure factor in k-space^i^"— . Exactly this situation is 
realized in the rare-earth magnet spin ice^. 

The purpose of this paper is to treat the quantum effects 
which arise from the tunnelling of the system from one dimer 



configuration to another, as described by Eq. ([T]). This tun- 
nelling introduces fluctuations in time of the gauge field A, 
which in turn give rise to an effective electric field 



E 



-dtA. 



(12) 



None the less, the total magnetic flux (j> = / dS • B through 
any plane in the lattice remains a conserved quantity, since the 
dynamics present in the quantum dimer model only permit a 
reversal in the sense of flux around a closed loop. 

This representation clearly has a lot in common with con- 
ventional electromagnetism and, following Ref.Hl we can use 
this analogy to write down a plausible long-wavelength action 
for the QDM on a diamond lattice 



S = J Sxdt [E^ - c^B^] 



(13) 



where c? > 0. This is the Maxwell action of conventional 
electromagnetism, and the system must posses linearly dis- 
persing "photons" (transverse excitations of the gauge field 
A) with speed of light c. These conditions were argued to 
be realized in a quantum U (l)-liquid state bordering the RK 
point in three-dimensional quantum dimer models on a bipar- 
tite lattice for /i < 1 ilt], leading to the conjectured phase 
diagram Fig.[T] An exactly parallel story can be told about the 
quantum loop model in three dimensionsi^. 

These "emergent" photons are the signature feature of the 
proposed J7(l)-liquid state, and offer a beautiful realization 
of Maxwell's laws in a condensed matter system. However, 
for the purposes of this study, we wish to emphasize that 
Eq. ( fTSI ) also contains information about the finite size scal- 
ing of the energy spectra. A flux through a cluster of vol- 
ume corresponds to an average magnetic field B — 
In the C/(l)-liquid state, this magnetic field is uniformly dis- 
tributed on the "coarse-grained" scale of the effective action 
Eq. ( [13] ). It then follows from Eq. (fTSl l that the energy differ- 
ence = E^ — Eq between the ground state of the zero-flux 
sector, and the lowest energy state of the sector with flux cj) 
scales as 



E. 



En = c'-^ . 



(14) 



In the thermodynamic limit and different flux sectors 

should be treated as degenerate ground states. A second finite- 
size prediction of Eq. ( fTsT i is that the ground state energy of 
the J7(l)-liquid state should have a finite size correction AE 
from the zero-point energy of photons which in leading order 
scales as 



AE{L) 
N 



A' 



dk 



A' 



dk 



k^ck 



c 

l4 



(15) 



where A' is a size-independent momentum cutoff. In this pa- 
per we make use of both Eq. ( fT4T ) and Eq. ( fTsl l to study finite 
size properties of Eq. ([T]i. 
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TABLE I: Cluster geometries used in simulations. The first column 
gives the short notation for cluster geometry used in the text; the 
second column, the number of diamond lattice bonds; the last three 
columns show the vectors gi which define the translation vectors of 
the cluster. The integer I defines the size of the cluster, within a given 
family. The [100] clusters have side of length of L = 2Z, and contain 
iV = 2L'-^ bonds. 





number of bonds 


gl 


g2 


g3 


[100] 


16/^ 


2; (1,0,0) 


21 (0,1,0) 


21 (0,0,1) 


[110] 


32/^ 


2Z (1,1,0) 


21 (0,1,1) 


2/(1,0,1) 


[111] 


64/^ 


2; (1,1,-1) 


2/(1,-1,1) 


2/ (-1,1,1) 



C. Clusters used in simulation 

Our evidence for the existence of a quantum /7(l)-liquid 
phase is taken from simulation of the quantum dimer model 
Eq. ^ on finite-size clusters of the diamond lattice with pe- 
riodic boundary conditions. For the size of a cluster use the 
notation in which N is the number of pyrochlore lattice sites, 
corresponding to the diamond lattice bonds. The periodic 
boundary conditions permit us to associate a set of (integer) 
topological quantum numbers (f> ~ {(f>i,(j)2,(t)3) with the flux 
through a set of orthogonal planes. This reduces the Hilbert 
space of the problem from the ?» 1.3^^^ different dimer con- 
figurations to a set of discrete sectors with given (f>. 

An important lesson from the exact diagonalization of small 
clusters is that the matrix elements of Eq. ^ do not connect 
all dimer configurations within a given value of <j>. None the 
less, the size of the largest blocks is still too large to permit the 
exact diagonalization of clusters with more that 128 diamond 
lattice bonds. These clusters are very small by the standards 
of the (implicitly) course-grained field theory Eq. (flST l, and 
quantum Monte Carlo simulations must therefore be used to 
determine the ground state of the model. 

The quantum Monte Carlo method we chose is GFMC 
(Green's function Monte Carlo), which constructs Monte 
Carlo updates using only the Hamiltonian dynamics of Eq. ([l]). 
This has the advantage that all quantum numbers (including 
magnetic flux) are preserved in simulation. However since not 
all of these quantum numbers are known a priori, great care 
must be taken to simulate in the appropriate subsectors of the 
Hilbert space. This is an issue which we discuss at length 
in Section |lVl Moreover, because an exponential number of 
dimer configurations contribute to the ground state wave func- 
tion of the /7(l)-liquid, quantum Monte Carlo simulations are 
themselves limited to clusters of up to 2000 diamond lattice 
bonds. 

In order to gain the most from finite size scaUng, we there- 
fore consider a range of different cluster geometries. We 
restrict ourselves to clusters which are compatible with the 
16-bond unit cell of the R-state, and group them into fami- 
lies of clusters with the same shape. These are summarized 
in Table |T] The main family of clusters used in our simu- 
lations consists of clusters with 2L'^ diamond lattice bonds 
(pyrochlore lattice sites) with side of length L = 4, 6, 16, 



which have the full (cubic) symmetry of the diamond lattice. 
We denote this family as [100], since the edges of the cluster 
are parallel to cubic lattice axes. The two other families of 
clusters considered are denoted [110] and [111]. All results 
are quoted for [100] clusters, unless specified otherwise. 



D. String defects and flux sectors 




FIG. 4: (Color online) Construction of the "string" defect (one mag- 
netic flux) configuration from the maximally flippable R-state con- 
figuration (assuming periodic boundary conditions): (a) two adjacent 
16-bond unit cells of the diamond lattice within the R-state, with flip- 
pable hexagons denoted by blue shading; (b) a pair of monomer ex- 
citations created by removing a dimer; (c), (d) and (e) monomers are 
separated. The dashed line shows the "string" of displaced dimers 
which connects them, (f) one of the monomers has traversed the 
cluster and the dimer is re-introduced. The resulting configuration 
has a single string excitation and a nonzero magnetic flux. 
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As discussed in Section IIIBI and Section III CI above, any 
given dimer configuration on the diamond lattice can be clas- 
sified according to its total "magnetic" flux 

<^ =(</)!, (/)2,</)3) (16) 

where each component of is a (topological) quantum num- 
ber conserved under the dynamics of the quantum dimer 
model Eq. ([T]i. 

In a finite size cluster, each component of flux (j>i serves 
as a winding number which takes on an extensive number of 
(discrete) values. These are bounded by the size of the sys- 
tem, and can be expressed as combinations of the four occu- 
pation numbers (ni, n2,n^, ^4) which count the number of 
dimers on each of the four fee sublattices of diamond-lattice 
bonds (equivalently, pyrochlore lattice sites). These occu- 
pation numbers are themselves invariant under the dynamics 
of the QDM. However the mapping between the occupation 
number representation (cf. Ref. [l2l) . and the flux represen- 
tation of these conserved quantities depends on the specific 
cluster geometry. 

For the (100) cluster series we obtain the following expres- 
sion: 

4>x = 2^*""^ ^ ^ ^ 

02 = ^("-1 + - "-3 - "4), (17) 

where the corresponding positions of the sites in the 4-site 
elementary unit cell of the pyrochlore lattice are numbered 

ri = (l,l,l)/4 

V2 = (-1,-1, l)/4 

rs = (-l,l,-l)/4 

r4 = (1,-1, -l)/4 

It follows that -L2/4 < 0, < L'^ /A for a cluster with N = 
2L^ bonds in the diamond lattice {i.e. the number of sites 
in the pyrochlore lattice). In the highly symmetric R-state 
(cf. Fig.|2]i ni = n2 = ?i3 = 714. This state therefore belongs 
to the zero-flux sector = (0, 0, 0). 

In order to test the predictions of the effective field the- 
ory, notably Eq. (O, we need to be able to construct other 
states in neighbouring flux sectors. The flux quantum numbers 
(01,02,9^3) are conserved under all local operations which 
connect configurations satisfying the dimer constraint. How- 
ever, we can construct states with finite magnetic flux, starting 
from the R-state, by creating non-local "string" defects which 
traverse the periodic boundaries of the cluster 

We do this following the procedure illustrated in Fig. HI 
When a dimer is removed from the lattice, two monomers 
(sites not touched by a dimer) are created. Subsequently, we 
can separate these monomers, creating a "string" of dimers 
which have been translated by exactly one bond. Moving a 
dimer by one bond reverses the sense of the magnetic field 
B associated with that dimer. We can therefore construct a 



state with a different flux by moving one of the monomers 
across the periodic boundary of the cluster in such a way that 
the two monomers meet and the dimer can be re-introduced. 
The resulting configuration fulfills again the zero-divergence 
condition on B, but the magnetic flux through the periodic 
boundary traversed by the monomer has changed by exactly 
one unit. For a cubic cluster this is the smallest magnetic 
flux possible in the system, and we choose convention where 
(/) = 1 in the new configuration. States with higher flux can be 
constructed by repeating the procedure above. 

The monomers created by removing a dimer are excitations 
with fractional quantum numbers which depend on the origi- 
nal microscopic model, e.g. they might be fractional charges 
e/2 ifTsll . or "spinons" carrying spin s = 1/2 ifisll . In field- 
theoretical terms, they are the magnetic monopoles of a U{\) 
gauge theory, and the central question addressed in this paper 
is whether the QDM on a diamond lattice can support a quan- 
tum U{\) liquid phase in which these monopoles are decon- 
fined. In a finite size system, deconfined monopoles must be 
free to traverse the periodic boundaries of the cluster, and so 
the string defects described above also provide a direct test of 
monopole confinement — a necessary condition for monopole 
deconfinement in a periodic cluster is that the ground states in 
flux sectors connected by strings should be degenerate [cf. 
Eq. (fT4l)1. We return to this idea below. 



E. Competing ordered phase 




FIG. 5: (Color online) Numbering convention in the definition of 
the order parameter vriR for the R-state (Eg.llSIl: (left) cubic cell of 
pyrochlore lattice with 16 sites, (right) corresponding diamond lattice 
cell with 16 bonds. 

In order to determine the phase diagram of the quantum 
dimer model on a diamond lattice (cf. Fig.[T]i, we need also to 
characterize the competing ordered "R-state" [Refs. [l^O- 
The R-state (illustrated in Fig. |2]i occupies the 8-site (16- 
bond) cubic unit cell of the diamond lattice (16-site on the py- 
rochlore lattice), and is 8-fold degenerate. The R-state breaks 
the inversion symmetry of the diamond lattice, its 8-fold de- 
generacy is best thought of as 2 x 4 states with opposite chi- 
rality. The 4 states with equal chirality can be related to each 
other with rotations or translations. In terms of the gauge field, 
this 8-fold degeneracy can be thought of as a 4-fold choice of 
111 axis (for a given tetrahedron), and a two fold choice of 
chirality about that axis. However for purposes of Simula- 
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TABLE II: Coefficients of bond occupation numbers within the six- 
dimensional irreducible representation used to define an order pa- 
rameter for the ordered R-state, following Eq. J18b . 
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lion it is easier to keep track of an order parameter defined 
in terms of the occupation numbers for dimers on diamond 
lattice bonds. 

We define this through a 6-dimensional irreducible repre- 
sentation of the diamond lattice point group as follows 



mR 



(18) 



where 



16 
«=1 



(19) 



Here mR,,, measures the occupancy of diamond lattice 
bonds/pyrochlore lattice sites, numbered as in Fig. |5l 
weighted by the factors given in Table HIl 

Naively, one might have expected the order parameter to 
have been given by the sum of projectors 

C=i 

where |R^) are the 8 possible ordered R-states 
(C = 1,2, ...,8). However, 



C=i 



N 

T 



measures the total number of dimers (equal to one quarter of 
the number of pyrochlore sites), and the only linearly inde- 
pendent combinations of |R(^) are those given in Table. Ull 



III. EVIDENCE FOR A QUANTUM [/( 1 ) -LIQUID PHASE 

A. Comments on simulations 

In this section of the paper, we calculate the ground state 
properties of the quantum dimer model Eq. ([T]) using the 
Green's Function Monte Carlo (GFMC)2^ method. This 
method has previously been applied very successfully to 



QDM's in two dimensions^^Ti^. We benchmark our GFMC 
simulations against results from exact diagonalization for 
small clusters, against a perturbation theory about the RK 
point for /^t — > 1 and, in AppendixiDl against expansions about 
a perfectly ordered R-state for /i < 0. We chose to use GFMC 
because it automatically conserves the flux quantum numbers 
which are central to our study of the QDM, and because it 
can provide numerically exact answers for the ground state 
(T = 0) properties of reasonably large clusters. 

GFMC works by calculating a set of running averages on a 
random walk through configuration space, using on Hamilto- 
nian matrix elements to move from one configuration to an- 
other. For large clusters and/or complicated problems, this 
random walk must be guided using a trial wave function, pre- 
viously optimized by a suitable variational calculation. In this 
sense GFMC can perhaps best be understood as a systematic 
way of improving upon a known variational wave function. 

In Appendix |C] we discuss a range of different variational 
wave functions which can be used to explore the ground state 
of the QDM on a diamond lattice. The most useful guide func- 
tion we have found to date is ; 

\'^'^Pi)4>.^ = exp[aA^f + /3tor + ^ lijnT;j] |^/'rk)0x (20) 

Here |V'rk)</).c is the equally weighted superposition of all 
dimer configurations within the maximally-connected subsec- 
tor with flux 4> ~ (0i, (/)2, 03 ), i.e. those which are connected 
to an R-state by string excitations and/or matrix elements of 
Eq. ([T]). The sum ij runs over pairs of hexagonal plaquettes 
within the diamond lattice, and r, = 0,1 is an Ising-Iike vari- 
able which takes on the value T,j = 1 when the hexagonal pla- 
quette i is "flippable". A'f ~ ^ - counts the total number 
of flippable plaquettes and mR is the order parameter associ- 
ated with the R-state, as defined in Eq. ( flST l. The variational 
parameters a, (5 and are minimized using a stochastic re- 
configuration algorithm^. For large lattices there are a many 
different pairs of hexagonal plaquettes. In practice we restrict 
ourselves to a set of up to 40 inequivalent 7^ . 

Our GFMC algorithm follows the prescription of [27 1 . Sim- 
ulations are typically performed using up to a few thousands 
of "walkers" in parallel. The transition probabilities between 
dimer configurations are calculated using the matrix elements 
from Eq. ([T]). If the guide wave function is exact (for example 
at the RK point, with all variational parameters set to zero) 
GFMC converges to the exact answer in a single step. Away 
from the RK point simulations typically converge after a few 
hundreds of reconfiguration steps. A "forward walker" tech- 
nique is used to calculate the average value of the order pa- 
rameter mR ifz^. 

For = 0, all simulations are started from a perfectly 
ordered R-state, and explore only dimer configurations con- 
nected with this. For more general 0, starting configura- 
tions are constructed using the prescription described in Sec- 
tion IIIDI In neither case is the GFMC technique ergodic, in 
the sense of exploring all dimer configurations which might, 
in principle, contribute to the ground state. However GFMC 
performed in this way does produce results which are repre- 
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sentative of the entire configuration space. This question is 
discussed at length in Section lTVl 



B. Death of an ordered state 

We begin by studying how the ordered R-state evolves as a 
function of /i. In Fig.|6]we present both VMC and GFMC sim- 
ulation results for the order parameter mp, for —0.2 < fJ, < I, 
in a series of cubic clusters of the form [100] which have 
the full symmetry of diamond lattice and have 2L^ bonds 
[cf. Table in. The linear dimension of the clusters ranges from 
L = 4 (128 bonds) to L = 10 (2000 bonds). 



0.4 




I , I . I , I . I , I . 

-0.2 0.2 0.4 0.6 0.8 1 

1^ 



FIG. 6: (Color online) Collapse of the ordered state: the order pa- 
rameter TTiR in the sector connected to the R-state. The results of 
VMC (dashed lines) and GFMC (solid lines) simulations are shown 
for series of clusters with cubic symmetry, and compared to exact 
diagonalization results for the smallest, 128 bond cluster. The inset 
shows finite size scaling of the order parameter at the RK point. 

For the smallest 128-bond cluster we find very close agree- 
ment between GFMC and VMC simulations, and perfect nu- 
merical agreement between GFMC results and exact diago- 
nalization calculations. This cluster is too small to be repre- 
sentative, but a point of inflection in ttir for /i « hints at a 
possible phase transition out of the R-state. For larger clus- 
ters, a first-order phase transition out of the R-state is clearly 
visible as a jump in the order parameter ?tir at a critical value 
of /i = /ic(i)- This jump is most pronounced in VMC cal- 
culations, where it occurs for smaller values of /i, i.e. deeper 
inside the ordered state. The jump in the ttjr shifts steadily to 
larger values of fi as the size of the system is increased. For the 
L ~ 6, 432-bond cluster, a jump in tor is observed in VMC 
for fic ~ 0.25, and in GFMC for /ic « 0.45. Meanwhile for 




FIG. 7: (Color online) Estimation of the critical value of /ic (oo) 
for the transition from the ordered R-state to the quantum U{1)- 
liquid, taken from GFMC simulation of three different families of 
clusters. Linear fits to fJ.c{L) as a function of 1/L* are made subject 
to the constraint that all data sets must converge to the same value 
/ic(oo) for l/L — > 0, regardless of the cluster geometry. We find 
/ic(oo) = 0.75 ±0.02. 



the largest L = 10, 2000-bond cluster, the GFMC simulations 
suggest that this transition takes place for fic « 0.7. 

For i-Lc{L) < /i < 1, niRliJ,) takes on a smaller, roughly 
constant value. However the finite size effects are large, with 
the value of tor decreasing as the system size is increased. 
Much larger systems can be simulated at the RK point /i = 1, 
where the form of the ground state wave function is known 
exactly. Various methods of simulating at the RK point are 
discussed in Section |IV] In the inset to Fig.|6l we show re- 
sults obtained using local updates within the configurations 
explored by GFMC simulations. We find that tor vanishes in 
the thermodynamic limit as TOR(/i = 1) ~ L~^/'^. 

Given the strong finite size effects visible in these data, it 
is important to make an estimate of iic{L oo), to ensure 
that the competing phase for ^ > fJ,c{L) does not collapse to 
a single point in the thermodynamic limit. The properties of 
the system for L ^ oo must be independent of the cluster ge- 
ometry, and so fic{L) should interpolate to the same value of 
/ic(oo) for all different families of clusters. In Fig. |2]we plot 
Hc{L) as a function of l/L"^ for clusters of the type [100], 
[110] and [111], as defined in Table U All data collapse to a 
single value fic{oo) = 0.75 ± 0.02. We note that the slightly 
different method of estimation used in Ref . [13 (fitting a sepa- 
rate line to each family of clusters) gave a very similar result 
Aic(oo) = 0.77 ±0.02. 

The fact that fic{oo) ~ lJ-c{L) ^ provides strong, al- 
beit indirect, evidence for the existence of the gapless photon 
excitations which are a signal feature of the quantum U{1)- 
liquid state. We can argue as follows : since the jump in 
Tor implies a first-order, zero-temperature phase transition, 
the finite size correction to /ic must follow from the finite size 
corrections to the ground state energy of the two competing 
phases. All excitations about the R-state are gapped, and finite 
size corrections to its ground state energy should therefore be 
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1^ 

FIG. 8: (Color online) Monopole string tension ki [Eq. 1211 for 
a range of ^ spanning the ordered R-state and proposed quantum 
(7(l)-liquid phase. Results are taken from VMC (dashed lines) 
and GFMC (solid lines) simulations for a series of clusters with cu- 
bic symmetry, and compared to exact diagonalization results for the 
smallest, 128 bond cluster. 



small. Assuming that the competing phase is a quantum U{1)- 
liquid with gapless photon excitations, the finite size correc- 
tions to its ground state energy will scale as ^ 1/ i^, as de- 
scribed by Eq. ( fTSl l. Provided that the ground state energies 
of both states evolve smoothly near fJ-dL) (which is known to 
be true from both exact diagonalization and quantum Monte 
Carlo studies), finite size corrections to /ic should also then 
scale as ^ l/L"*. 

On the strength of the simulation results for ttir presented 
above, it seems very reasonable to conclude that the phase 
diagram proposed in Fig.[T]is correct, with a single first-order 
phase transition from an ordered R-state to a quantum U(l)- 
liquid occurring for Hc ~ 0.75 ± 0.02. In what follows, we 
present a variety of other evidence that this is indeed the case. 



C. Collapse of monopole string tension 

If the state occurring for fJ-dL) < fi < 1 is indeed a 
quantum C/(l)-liquid, it should possess gapped, deconfined 
monopole excitations'*. We can study the motion of these 
monopoles indirectly through the thought experiment used 
to construct configurations in different flux sectors in Sec- 
tion IIIDI As illustrated in Fig. HI separating the two 
monomers (monopoles) created by removing a dimer, de- 
stroys flippable hexagons along the "string" of dimers con- 
necting the monomers. In the crystalline state, where the po- 
tential energy dominates, one can expect that there is a finite 
energy cost of such an operation, proportional to the length of 
the string. Where the string defect threads a periodic bound- 
ary of the cluster, its energy cost can be calculated as the gap 
to the ground state with one additional unit of flux. 

We therefore define the "string tension" as the energy dif- 
ference between the ground state in the zero-flux sector (Eq) 
and a sector with a single string defect (Ei), divided by the 



length of the string (L) 



For the string defects considered, L is simply the linear di- 
mension of the cluster This string tension is an indirect mea- 
sure of the attractive, confining force between monopole ex- 
citations. It must therefore vanish in the deconfined quantum 
C/(l)-liquid, but will be finite in the any ordered state. 

In Fig. [8] we plot simulation result for the monopole string 
tension ki as a function of /i for the same set of clusters and 
range of parameters as Fig. |6] Deep inside the ordered state, 
for negative values of i^l, classical effects dominate. A single 
string defect in a perfectly ordered R-state destroys 2L flip- 
pable hexagons at a potential energy cost of 2L|/i|. This lin- 
ear relationship between Ki and for < is clearly visible 
in GFMC simulation results for the larger clusters. However 
the string tension drops rapidly as we approach the transition 
from the R-state into the competing phase at fj, ^ f^ciL), and 
for large clusters it is essentially zero for ndL) < /i < 1. 

The coincidence between the jump in the order parameter 
rriR [Fig.|6l and the collapse in the string tension ki [Fig.[8l 
lends yet more support to the proposed form of the phase dia- 
gram for the QDM on a diamond lattice [Fig.[T]. In particular, 
the absence of a string tension for fJ,c{L) < /i < 1 rules out 
the possibility that the competing phase is a different ordered 
state, for example a three-dimensional analogue of the res- 
onating plaquette phase found in quantum dimer and quantum 
ice models in two dimensions^i^2i^. However we still need to 
rule out the possibility that the competing phase is a different 
form of quantum liquid, by confirming explicitly that it is in- 
deed the quantum C/(l)-liquid we have been seeking. This is 
accomplished below. 

D. Birthof a (7(1) -liquid 

Within the field-theoretical scenario for a three dimensional 
QDM on a bipartite lattice, the quantum ?7(l)-liquid phase 
"grows" out of the RK point^ '^. Exactly at the RK point, the 
ground states in all flux sectors are degenerate. Away from the 
RK point, in a finite size system, the ground state in the zero- 
flux sector has the lowest energy, and the (finite-size) energy 
gaps to the other low-lying flux sectors scale as ^ c^clP /L 
[cf. Eq. (fT4]i1. By construction, c vanishes at the RK point. 

Since the exact ground state is known at the RK point, we 
can calculate the finite size gaps (and implicitly the speed 
of light c) within perturbation theory about the RK point. 
Writing 

■Hqdm = ^RK - S^iAff (22) 

where Hrk is defined by Eq. (|5]l, Sfi = 1 — /i, and the operator 
Aff is defined by Eq. (O, we find 

A0,MC = AN^^McSfi + 0(6(1^) (23) 

where 

AN^MC = {M)1%c - m^MC (24) 
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FIG. 9: (Color online) Reduction in the average number of flippable 
plaquettes AN^, as calculated in simulations at the RK point, (a) 
—ANtj, plotted as a function of for cubic clusters of linear dimen- 
sion L = 4 (432 bonds) to L = 16 (8192 bonds), (b) The same 
data set plotted as —AN^/(j>'^ vs to extract the leading depen- 
dence on (j>. (c) Finite-size scaling of the the coefficient of (/>^ found 
from fits to —AN^ for clusters of linear dimension L [cf. intercept 
to ordinate axis in (b), above]. The linear collapse of the coefficient 
of (p^ with 1/L is consistent with the prediction of the effective field 
theory for a quantum U{1) liquid Eq. l ll4b . 



is the difference in the average number of flippable plaque- 
ttes between the zero-flux sector and the sector with flux (f>, 
calculated at the RK point, within the maximally-connected 
subsectorof the Hilbert space, Ac = MC. The maximally flip- 
pable R-state belongs to the zero-flux sector, and we therefore 
anticipate that AA^^ < 0. If, as supposed, a quantum U{1)- 
liquid grows adiabatically from the RK point, then we must 
find 



AN,. 



L 



(25) 



In Fig.|9]we show results for AA^^ obtained from simula- 
tion of a series of [100] cubic clusters with linear dimension 
ranging from L — 6 (432 bonds) to L = 16 (8192 bonds), for 
flux (f) = (0,0,0), plotted as a function of (fP. Simulations 
were performed within the maximally-connected subsector of 
states connected to the R-state for = 0, and within the sub- 
sector associated with a given number of string defects for 
higher values of flux. Due to large finite size effects, the ex- 
pected linear relationship between —AN^ and 0^ is only vis- 
ible for clusters with more than 2000 bonds. However the co- 
efficient of 0^ can be extracted from a polynomial fit to AA^^ 




(b) 



I 



0.1 



0.05 



1 1 




1 1 










♦ 


perturbation theory 








H = 0.99 






• 


H = 0.95 






• 


[1 = 0.9 




I 1 




1 1 



10 



20 



30 



40 



50 



FIG. 10: (Color online) Comparison of the energy gaps found 
in GFMC simulations with the predictions of a perturbation theory 
about the RK point [Eq. )22b ]. for a range of fi within the pro- 
posed quantum [/(l)-liquid phase, (a) Energy gap normalized 
to —Sfi — fi — 1, to extract the leading dependence on /x, and plot- 
ted as a function of 0^ [cf. Eq. (I14H. (b) Energy gap A^ normalized 
to —5/1 X ff)^, to extract the leading dependence on both fi and (j>, 
plotted as a function of (p^ [cf. Eq. J14b]. All results are for a 1024- 
bond cluster with cubic symmetry. 



as a function of 0, at given L. This is found to be proportional 
to 1/L, as shown in the inset to Fig. |9l From these results we 
obtain 



0.6 X 5/i 



(26) 



(in the present units) confirming that an incipient quantum 
U (l)-Iiquid is present at the RK point. This result is specific 
to the QDM on a diamond lattice. We note that simulations of 
the QDM on a cubic lattice yield quite different results^. 

In Fig.[TO]we compare the results of the perturbation theory 
Eq. ( I23I 1 with GFMC calculations made for a set of parameters 
bordering on the RK point, for a cluster of 1024 diamond lat- 
tice bonds. The energy gaps obtained from GFMC are plotted 
as Aif,/{1 — fi) so to compare with perturbation theory, and 
as a function of cf)^ so as to extract the leading behaviour of 
the t/(l)-liquid phase. Both sets of data are found to show the 
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FIG. 1 1 : (Color online) First order quantum phase transition from 
the ordered R-state into the quantum [/(l)-liquid, as imaged in the 
finite size spectrum of a 1024-bond cluster GFMC simulation results 
for the energy gap are plotted as Arf,/[(j>'^{l — fi)] so as to collapse 
the spectra for the ?7(l)-liquid phase. The phase transition is visible 
as an abrupt qualitative change in the form of the spectra for /j, — 0.7. 
The results of perturbation theory about the RK point (black line and 
points) are shown for comparison. 



same (fy^ dependence on flux sector, confirming the existence 
of a quantum J7(l)-liquid bordering the RK point. 



E. A phase transition in a spectrum 

We are now in a position to draw a set of strong conclusions 
about the zero-temperature phase diagram of the QDM on a 
diamond lattice (cf. Fig.[l]i- For 

At < /ic = 0.75 ± 0.02 

the ground state of the model is an ordered R-state. For 
/ic < M < 1 it is a quantum [/(I) -liquid with linearly dispers- 
ing photon excitations and deconfined magnetic monopoles. 
The quantum phase transition between these two phases at 
/i = /ic is first order. For /t — > 1_, the quantum J7(l)-liquid 
is adiabatically connected with the RK point. 

In Fig.[TT]we present the results of GFMC simulations of a 
cubic 1024-bond cluster, for a range of parameters < /t < 1 
spanning both the R-state and the quantum L/(1)-Iiquid. We 
plot the energy gaps A^, rescaled as Acf,/[(fP{l — fi)] so as to 
collapse all data within the t/(l)-liquid phase. Deep within 
the ordered phase the should be proportional to both the 
linear dimension of the cluster and the flux (number of string 
defects), i.e. 



A 



L X 



A plot of A^ / 0^ at fixed L should therefore show a clear 
distinction between a t/(l)-Iiquid (A^/c/)^ ~ const.), and an 
ordered phase (A^/^^ ^ 1/0)- For this cluster size a clear 
division is observed between /t > 0.7 (liquid) and /t < 0.6 
(ordered). This is unambiguous evidence of a phase transition 
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FIG. 12: (Color online) (a) Block-structure of the Hamiltonian 
for the quantum dimer model on a diamond lattice, for an 128- 
bond cluster with a fixed value of the flux (ji = (0, 0, 0), following 
Table mil Connected dimer configurations break up into different 
sub-ensembles. These can be divided into two groups, according 
to whether their transition graphs contain an odd or even number 
of loops. Both even and odd subsectors contain isolated configu- 
rations which do not participate in the dynamics of the model, and 
do not contribute to the ground state for /t < 1. (b) Contrasting 
block-structure of the Hamiltonian for the quantum dimer model on 
a square lattice for a finite size cluster with a fixed value of the flux 
= (0, 0). All dimer configurations which contribute to the ground 
state for /t < 1 belong to a single connected block. The remaining, 
isolated, dimer configurations do not participate in the dynamics of 
the model. 



from a linearly confining phase to a quantum [/(l)-liquid, im- 
aged in a spectrum. The abrupt, qualitative change in the spec- 
tra for /i = 0.6 confirms that this phase transition is first-order 
in character. 



IV. ERGODICITY AND CONSERVATION LAWS WITHIN 
A GIVEN FLUX SECTOR 

A. Lessons from exact diagonalization 

Since GFMC uses on matrix elements of the QDM Hamil- 
tonian [Eq. ([T])] to connect different dimer configurations, it 
is guaranteed to respect the flux quantum numbers defined 
in Section HID I However, not all dimer configurations with 
the same magnetic flux are connected by matrix elements of 
the Hamiltonian. Clearly, therefore, GFMC simulations for a 
given flux are not ergodic, in the sense that they do not explore 
all dimer configurations which might, in principle, contribute 
to the ground state for that flux sector. It is therefore very 
important to consider how the Hilbert space for a given flux 
(f) breaks up into different subsectors, and to ensure that sim- 
ulations are carried out in the appropriate sub-sector(s). To 
this end, much can be learnt from the exact diagonalization of 
small clusters. 

In total there are 115, 150,848dimer coverings ofthe L = 2, 
128-bond, cubic cluster. Of these, 11,835,352 have zero flux 
{(j) = 0). These zero-flux configurations include 1440 isolated 
configurations with no flippable plaquettes, and 3072 config- 
urations with a single flippable plaquette which connect only 
to a single other configuration. The remaining 11,830,840 
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configurations can be grouped into a small number of subsec- 
tors which are connected by the matrix elements of Eq. ([T]). 
These are listed in Table |III1 and illustrated schematically in 
Fig.[T2ja). This situation should be contrasted with the much 
simpler structure of the Hilbert space in the QDM on a square 
lattice, where all connected dimer configurations for a given 
value of flux = {(px, <t>y) are associated with a single block 
of the Hamiltonian — cf. Fig.[T2lb). 




TABLE III: Subsectors of the Hilbert space of the QDM on a dia- 
mond lattice with zero flux {cj> — 0) for a cubic, 128-bond cluster. 
Subsectors are listed in order of increasing dimension, with isolated 
configurations and subsectors consisting of only two configurations 
omitted. The largest, maximally-connected, subsector (A) contains 
the R-states. Further subsectors can be obtained through cyclic per- 
mutations of dimers on closed loops — the 8-link loop operation 
shown in Fig.|13|provides access to subsector B. Subsectors can be 
classified according to the number and type of loops in their transi- 
tion graph. 



subsector 
label 

Ac 


multiplicity 


dimension 


number of 
loops in 
transition 
graph (mod 2) 


longest 
irreducible 
loop in 
transition graph 


A 


1 


7794744 





2 


B 


1 


3468032 


1 


8 


C 


3 


131584 





10 


D 


3 


45184 





14 


E 


8 


3376 





14 


F 


192 


40 
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FIG. 13: (Color online) 8-link "loop" used to connect dimer config- 
urations within subsectors A and B of the Hilbert space of the QDM 
on a diamond lattice [cf. Table IIIIL (a) 8-link loop of bonds within 
the cubic unit cell of the diamond lattice. Three distinct 8-link loops 
of this type are possible, each of which appears octagonal when pro- 
jected into an appropriate [100] plane. Where exactly 4 bonds are 
occupied by dimers, this loop is "flippable", i.e. a new state obey- 
ing the dimer constraint can be obtained by cyclicly permuting those 
dimers around the loop, (b) equivalent 8-link loop of sites within the 
cubic unit cell of the pyrochlore lattice. 




(a) 





For this 128-bond cluster, about two-thirds of the configura- 
tions with = belong to a maximally-connected subsector, 
labelled A in Table |llll which includes the R-states. The vast 
majority of the remaining configurations — about one-third 
of the total — are found in one other subsector, labelled B in 
Table HU] Subsector B can be accessed from the maximally 
connected subsector by acting on a randomly chosen configu- 
ration with the update shown in Fig. [13] — the cyclic permuta- 
tion of dimers around a closed, 8-link loop. Interestingly, this 
8-link permutation is the sub-leading term for dimer dynamics 
in the derivation of the QDM from Eq. (|7| or Eq. dS]). 

Subsectors C, D, E and F cannot be accessed from sub- 
sector A using the 8-link loop alone. However they can 
be reached by cyclically permuting dimers within suitably 
chosen configurations belonging to subsector A, on loops of 
length 10, 14, 14 and 12, respectively. Since the number and 
length of possible such loops grows with the size of the sys- 
tem, we anticipate that the number of subsectors in the Hilbert 
space grows with system size. We put these observations on a 
more formal footing below. 



FIG. 14: (Color online) Schematic representation of all possible 
changes in the topology of a transition graph which follow from the 
cyclic permutation of dimers on a hexagonal 6-link loop. Red lines 
denote loop segment with a single dimer. Blue lines denote a loop 
segment with an odd number dimers. All loop segments, whether red 
or blue, start and finish on sites belonging to different sublatticies 
of the (bipartite) diamond lattice, denoted here by filled or empty 
circles. The number of loops associated with each section of the 
transition graph is shown in the middle of the hexagonal plaquette. 
This number changes by either or 2. 



B. A hidden quantum number 

The way in which the Hilbert space of the QDM on a dia- 
mond lattice breaks up into distinct blocks suggests the exis- 
tence of further, conserved quantities besides magnetic flux cf). 
We can gain more insight into this problem, and understand 
the action of the 8-link loop shown in Fig. dTsl ) by studying 
transition graphs. 
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FIG. 15: (Color online) Schematic representation of all possible 
changes in the topology of a transition graph which follow from the 
cyclic permutation of dimers on an 8-link loop. The number of loops 
of the transition graphs changes by either 1 or 3. 



A transition graph is the set of directed loops formed by 
overlaying two dimer configurations on a bipartite lattice, and 
drawing arrows from A to B sublattice on those bonds which 
are occupied by a dimer in the first configuration, and from B 
to A sublattice on those bonds which are occupied by dimers 
in the second configuration"*. Cyclicly permuting the dimers 
on the loops of the transition graph converts the first configu- 
ration into the second (and vice versa). By construction, these 
loops cannot intersect one another — although one loop may 
be threaded through another in a three-dimensional lattice. 

A pair of dimer configurations are connected if (and only 
if) we can undo all the long loops in their transition graph us- 
ing only matrix elements of Eq. ([l]) — cyclic permutations 
of dimers on a 6-link hexagonal plaquette — so that only 
trivial "loops" occupying a single bond remain. These tiiv- 
ial loops have length 2, and the maximum number of trivial 
loops achievable is equal to the number of dimers in the sys- 
tem — an even number for all of the clusters considered in this 
paper In Fig. [14] we enumerate the three possible, topologi- 
cally distinct, changes in the transition graph associated with 
the cyclic permutation of dimers on an hexagonal plaquette. 
In each case the total number of loops in the transition graph 
changes by an even number — either or 2. 

Cyclic permutations of dimers on an 8-link loop have a 
somewhat richer structure, illustrated in Fig. [Ts] In this case, 
the total number of loops in the transition graph changes by 



an odd number — either 1 or 3. Consequently, the action of 
8-link loop cannot be undone using matrix elements of the 
Hamiltonian [Eq. ([B] alone. 

More formally, we can associate a unique transition graph 
with each distinct dimer configuration by choosing one partic- 
ular dimer configuration as reference (for present purposes , an 
R-state is a convenient choice). The number of loops in this 
transition graph, modulo 2, is then a conserved quantity within 
the dynamics of the QDM Eq. ([T]). Furthermore, it must be 
even in the subsector connected with the R-state. In Table Hill 
we classify the all of the different subsectors of the zero-flux 
sector of the 128-bond cubic cluster containing more than two 
configurations according to this new, Ising, quantum number. 

Examination of transition graphs also provides some insight 
into the other, smaller, subsectors of the Hilbert space, la- 
belled C, D, E and F in Table Hill Empirically, we find that 
these can be accessed from the maximally-connected subsec- 
tor containing the R-states (subsector A) through the cyclic 
permutation of dimers on suitably chosen loops of length 10, 
12, and 14. Closed loops in the transition graph can generally 
be contracted through the repeated permutation of dimers on 
6-link loops of the type found in the Hamiltonian Eq. ([T|). By 
construction, all non-trivial loops in the transition graph for 
two configurations within subsector A can be removed in this 
way, leaving only trivial loops of length 2. The different sub- 
sectors B — F listed in Table|lll|can therefore be classified ac- 
cording to the longest irreducible loop in their transition graph 
with a (suitably chosen) state within subsector A. The lengths 
of the irreducible loops are listed in Table Hill 

It is interesting to note that the Ising quantum number asso- 
ciated with the number of loops in the transition graph is even 
where the longest irreducible loop is of length 2, 10 or 14, 
and odd where the longest irreducible loop is of length 8 or 
12. We infer (but have not proved) that topological selection 
rules of the form illustrated in Fig. [T4| apply for updates based 
on loops of length 6, 10, 14. . . , and of the form illustrated in 
Fig.fTSlfor updates based on loops of length 8, 12, 16. . . Cyclic 
permutations on loops of length 6, 10, 14. . . are precisely the 
form of dynamics permitted for (spinless) fermionic models 
of the form Eq. while bosonic models also permit cyclic 
permutations on loops of length 8, 12, 16. . . We therefore 
speculate that the number of loops in the transition graph will 
be a conserved quantity (mod 2) for the most general QDM 
Hamiltonian derived from a fermionic model on a pyrochlore 
lattice, but not for the equivalent bosonic or spin model. 

Clearly, more questions could be asked about the non- 
ergodicty of the QDM on a diamond lattice. What, for ex- 
ample, are the the conserved quantities which distinguish the 
smaller subsectors of the Hilbert space ? And how do the 
number, multiplicity, and dimension of these subsectors grow 
with the size of the system ? For the purposes of this paper, 
however, the most pressing need is to understand the impact of 
this non-ergodicity on the simulations described in SectionHU 
It is this to which we now turn. 
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C. Can we trust simulations ? 

It is clear from the analysis above that QDM on a diamond 
lattice is not ergodic within a given flux sector The analy- 
sis which leads to the phase diagram Fig [T| rests on simula- 
tions performed within the (maximally connected) subsector 
of states connected to the R-state by the matrix elements of 
Eq. ([TJ, and on the systematic construction of configurations 
with finite flux using "string" defects' [cf. Fig. ©I. Before 
we can have full confidence in these results we must therefore 
address the following question ; 

What bias do we introduce into our results by simulating 
only within these subsectors of the Hilbert space ? 

The enumeration of states for an 128-bond cluster confirms 
the common-sense expectation that we can access the largest 
subsector of the Hilbert space for zero flux by seeding simula- 
tions from an R-state. Moreover, we now understand how to 
extend these simulations to configurations with an odd num- 
ber of loops in their transition graph, through use of an 8 -link 
loop. However there remain a non-vanishing fraction of dimer 
configurations — about 5% for the 128-bond cubic cluster 
[Table Unl — which can only be reached by other, more com- 
plex updates. 

It seems reasonable to expect that these unsociable con- 
figurations will (a) contribute little to the ground state wave 
function and (b) constitute a vanishing fraction of the total 
number of dimer configurations in the thermodynamic limit. 
However this reasonable expectation must be approached with 
some caution, not least because, in the absence of a full un- 
derstanding of conserved quantities, it is hard to assess how 
these neglected corners of the Hilbert space grow with system 
size. 

We can most easily get a handle on the consequences of 
non-ergodicity at the RK point, where simulations for ground 
state properties amount to performing classical averages over 
all dimer configurations within a given flux sector These can 
be performed in several different ways 



1. As in Section ITlIDI using "local" updates based on the 
matrix elements of the Hamiltonian. 

2. Using a "loop" algorithm developed for the simulation 
of ice model&i^i^. 

3. Using a "worm" algorithm based on movement of 
monomers^ 

In the "loop" algorithm, a directed random walk along 
bonds is used to construct a self-intersecting path in which 
all flux arrows associated with dimers have the same sense. 
The dimers which lie along this path are each cycled one link 
further along, reversing the sense of the flux around the closed 
"loop", while preserving the dimer constraint. 

In the "worm" algorithm, a single dimer is removed, and the 
two monopoles created are allowed to diffuse freely around 
the lattice until they meet, at which point the missing dimer 
is re-introduced. We have checked explicitly that "loop" and 
"worm" simulations lead to identical results for the QDM on 
the diamond lattice, and concentrate here on comparing re- 
sults obtained using "local" and "loop" updates. 
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FIG. 16: (Color online) Histogram of the number of flippable pla- 
quettes in different sets of dimer configurations, calculated at RK 
point, (a) Results for an 128-bond cubic cluster, calculated within 
in the entire zero-flux sector (blue, shaded columns), and in the sub- 
sector connected to the R-state (unshaded columns), (b) Equivalent 
results a for an 1024-bond cluster 



A necessary condition for the existence of a quantum U{1)- 
liquid in a QDM is that the average number of flippable 
plaquettes in the dimer configurations belonging to a given 
flux sector scale as AN^ = Nf(<j)) - iVf(O) ~ ~(t>'^/L 
[Eq. (l25Tll. In Fig. \Wia) we compare the probability distri- 
bution for the number of flippable plaquettes iVf calculated 
for all dimer configurations within the zero-flux sector (green, 
shaded columns) with that calculated within the maximally 
connected subsector containing the R-state (columns without 
shading). Identical results were obtained both within simula- 
tion, using "loop" updates to explore all dimer configurations, 
and "local" updates to access to the maximally-connected sub- 
sector, and within exact diagonalization, by explicit enumer- 
ation of all dimer configurations. The probability distribution 
for the number of flippable plaquettes is markedly different in 
the two cases, and the mean value of A^f clearly higher for the 
subset of configurations connected with the R-state. 

At first sight this disagreement might seem rather alarming, 
but it is simply a finite-size effect. In Fig. [T6lb) we present 
the results of an identical study for a larger, 1024-bond clus- 
ter The same, Gaussian, probability distribution with mean 
(Nf) = 73.27 and variance aN, = 6.21 is found for both 
the full set of dimer configurations sampled by the "loop" al- 
gorithm, and the maximally-connected subsector explored by 
the "local" updates. Similar gaussian probability distributions 
are found in sectors of the Hilbert space with finite flux, and 
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FIG. 17: (Color online) /i-dependence of the order parameter ttir 
for the R-state, within two different subsectors of the Hilbert space: 
starting from the R-state in flux (j) = subsector (open symbols) 
and in the subsector that is connected to it by an 8-site loop (solid 
symbols). Results are calculated using Variational Monte Carlo sim- 
ulation for cubic clusters with 432 (triangles) and 1024 (circles) 
diamond-lattice bonds. 

are again independent of the method of simulation. 

There are two obvious possible explanations as to why the 
same probability distribution should be found for averages 
calculated at the RK point in the full set of dimer configu- 
rations, and those calculated within its maximally-connected 
subsector One is that the relative size of the maximally- 
connected subsector grows with system size, so that it comes 
to dominate all averages in the thermodynamic limit. The 
other is that all large blocks of the Hilbert space (within a 
given flux sector) are statistically very similar This might oc- 
cur if, for example, configurations in different subsectors dif- 
fered only by some local defect which could not be removed 
using matrix elements of the Hamiltonian, but which had no 
impact on the bulk properties of the state. 

To distinguish between these alternatives we can exam- 
ine how quantum Monte Carlo simulation results in the two 
largest subsectors of the Hilbert space depend on the subsec- 
tor Ac in which they are carried out. In Fig. [17] we present 
the results of a Variational Monte Carlo study of the order 
parameter {mR)\^, for a cubic 432 and 1024-bond clusters. 
Simulations were carried out within 

i) the subsector with flux (p = containing the R-state; 

ii) the subsector with flux </> = obtained by acting with an 
8-link loop update on a randomly chosen state from the 
subsector containing the R-state. 

Within the crystalline, ordered phase, an 8-loop update re- 
duces the TOR relative to the subsector containing the R-state. 
This reduction in order parameter scales as — as would 
be expected for a local defect into an ordered state — and van- 
ishes in the thermodynamic limit. Meanwhile the energy cost 
of the 8-loop update is a constant in the ordered state. This 
behaviour should be compared with the topological string de- 
fect used to change flux sector acting in the R-state ordered 



FIG. 18: (Color online) The difference in the average number of 
flippable plaquettes ATV^ [Eq. (I24H. calculated at the RK point for 
dimer configurations in different subsectors of the Hilbert space. 
Green circles denote results for configurations connected to the R- 
state by matrix elements of Eq. ([TJ and string defects. Maroon dia- 
monds denote results for configurations containing an additional 8- 
link "loop" defect. All results are for a cubic 1024-bond cluster. 



phase, whose effect on the order parameter scales as 1/ L^, 
and which creates a finite size energy gap oc L. 

In marked contrast, the simulation results for (mR(/i)), in 
the quantum U (l)-liquid phase are independent of the subsec- 
tor within which they are calculated. Intuitively, we imagine 
that as the R-state order melts into the t/(l)-liquid phase, so 
defects in the order parameter melt, too, leaving only an invis- 
ible partition in the Hilbert space. 

We have repeated this analysis for more general local de- 
fects in the R-state, seeding simulations from configurations 
generated using randomly chosen loops with 10, 12, 14, 16, 
18 and 20 links, applied states within the subsector contain- 
ing the R-state. After simulation starting from a configuration 
obtained using one 10-, 14- or 18-link loop we recover the 
same values of (mR(/i)) as in the maximally connected sector 
containing the R-state. Meanwhile simulates started from a 
configuration obtained using an 12-, 16- or 20-link loop ex- 
hibit the same values of (mR(/i)) as those in the second large 
subsector accessed using the 8-link loop. 

This is consistent with our analysis of the transition graphs: 
the 8-loop cannot be undone by Hamiltonian matrix updates 
and so we remain in a subsector not connected to R-state, 
for which the average value of the order parameter has to be 
lower. Loops of length 10, 14, 18... can typically be reduced 
to 6-Iink loops so that the resulting configuration remains in 
the maximally connected sector, while loops of length 12, 16, 
20... can be reduced to an 8-link and 6-links loops, leading to 
the second large subsector Although we do not have a pre- 
cise understanding of the origin of the additional subsectors 
found in the 128-bond cluster [Table HiH . from our simulation 
for larger systems we tentatively conclude that the statistically 
accessible configurations with an even(odd) number of loops 
in their transition graph are of the type A and B (i.e. connected 
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to the R state directly or by a single 8-link loop). 

It is also interesting to compare the changes in the liquid 
phase ground state energy in detail. The ground state energy 
in all subsectors must be zero at the RK point. In the liquid 
phase bordering on the RK point, the gaps between states with 
different magnetic flux should scale as ^ clP / L [Eq. (fT4b 1 . 
This in turn implies that, at the RK point, AiV^ ~ —cfr' /L 
[Eq. dZSTll. In Fig.[T8]we present the results of simulations car- 
ried out in maximally-connected subsector for each 0, and in 
the subsector connected to it by a single, local, 8-link "loop" 
defect. The differences in the values of AiV^ obtained in the 
two subsectors are smaller than the error bars on the simula- 
tion. 

This analysis represents only a partial solution to the prob- 
lem of how to understand different subsectors within the 
Hilbert space of a three-dimensional QDM. However for the 
purposes of this paper, we gain confidence that GFMC simu- 
lations performed in the maximally-connected subsector, fol- 
lowing the prescriptions of Section and Section HIH do give 
reliable estimates of the ground state properties of the QDM 
on a diamond lattice. 

We can trust these estimates in the ordered, R-state, be- 
cause the alternative subsectors correspond to excited states, 
and not to the ground state. And we can trust the estimates in 
the quantum U (l)-liquid phase because in that case all statis- 
tically significant subsectors contribute equally to the ground 
state. Thus, while the QDM on a diamond lattice is not er- 
godic within a given flux sector, this lack of ergodicity has no 
pathological consequences for the ground state phase diagram 
as we have calculated it. Our initial analysis suggests that the 
same may not be true of the QDM on a cubic lattice^. 



V. CONCLUSIONS 

In this paper we set out to show that a three-dimensional 
quantum dimer model (QDM) — the QDM on a bipartite dia- 
mond lattice — could support a quantum liquid ground state, 
whose excitations include linearly dispersing photons and de- 
confined magnetic monopoles of an underlying U{1) gauge 
theory. We have identified the order parameter associated 
with the competing, ordered "R-state" and, through use of 
quantum Monte Carlo simulation, have shown that it is finite 
only where the ratio ji of kinetic to potential energy terms in 
the QDM Hamiltonian is less than /i^ = 0.75 ± 0.02. We 
have further shown that (a) the string-tension associated with 
separating a pair of monomers (magnetic monopoles) van- 
ishes for i^L > and (b) the finite-size energy spectrum for 
/ic < /^t < 1 is consistent with the predictions of the Maxwell 
action of conventional quantum electromagnetism. 

We conclude that the QDM on a diamond lattice undergoes 
a first-order phase transition from an ordered R-state into a 
quantum U{\) liquid phase at /i = /ic. This liquid phase 
occupies an extended region of parameter space, terminating 
at the RK point i^l ~ \. This completes the analysis outlined 
in our earlier Letter [Ref. [JU, and confirms the validity of the 
phase diagram Fig.[Tl previously proposed in Refs l8lfl2l on the 
basis of field-theoretical arguments. 



We have also explored how the Hilbert space of the QDM 
on a diamond lattice breaks up into different, disconnected, 
subsectors, and identified a new, Ising, quantum number as- 
sociated with the number of loops in the transition graph of a 
given dimer configuration. We argue that the lack of ergod- 
icity within any given flux sector does not have pathological 
consequences for the ground state phase diagram, as found 
by quantum Monte Carlo simulations. We further show that 
the same ground state phase diagram is obtained, regardless 
of whether the dimers are treated as bosonic or fermionic ob- 
jects. 

Taken together, these results would seem to set the exis- 
tence of a quantum U{1) liquid phase in the QDM on a di- 
amond lattice beyond reasonable doubt. However the very 
fact that such a state exists in this model begs all manner of 
other questions. Do other quantum dimer models on bipar- 
tite lattices in three dimension — for example the cubic lat- 
tice — also support quantum U{\) liquids, as conjectured in 
Ref. H ? What are the nature and consequences of the ad- 
ditional quantum numbers associated with the non-ergodicity 
of three-dimensional quantum dimer models ? Which other 
models can support a quantum U{1) ground state, and for 
what range of parameters ? An obvious candidate in this case 
is the quantum loop model studied in Ref. Here we find 
evidence for a quantum U{1) liquid for —0.5 < /i < 1, which 
will be presented elsewhere^. 

The greatest excitement, however, would attend to finding 
a quantum U{1) liquid phase in experiment, and identifying 
its deconfined, fractional excitations. Here more quantita- 
tive analysis is needed, both to pin down the characteristics 
of these excitations, and to bridge the gap between idealized 
models such as Eq. ([T]i, and real materials or artificial assem- 
blies of cold atoms. Since the coherence temperature of the 
quantum U{1) liquid will be low in either case, both finite 
temperature effect&i^, and the consequences of competing in- 
teractions, need to be explored in more detail. 
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Appendix A: Absence of a fermionic sign problem 

All results for Hqom are obtained under the assumption that 
the dimers are bosonic objects, i.e., exchanging two dimers 
does not yield a phase. However, as we show in this Appendix, 
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the same effective model can be used to describe the fermionic 
case. 




FIG. 19: (Color online) Labeling of the sites of the pyrochlore lattice 
by enumerating one of the two types of tetrahedra and four sublat- 
tices on each tetrahedron. This labeling allows us to show that the 
Hamiltonian Hqdm has no fermionic sign problem even if the under- 
lying particles are fermions. 

Let us now assume that the dimers are objects which are 
created (annihilated) by fermionic operators ^ (c^ ^). The 
index i E {1, . . . , n} refers to a tetrahedron and the index a G 
{1, 2, 3, 4} refers to the site on that tetrahedron (see Fig. 
Since the diagonal part Afj of the Hamiltonian in -ffgoM is 
not affected by changing to fermions, we only focus on the 
ring-exchange processes. These can be written in terms of the 
fermionic operators as 

(AO 

where the sum is taken over all hexagons, and ia < h < ic- 
We thus choose only one of the twelve different ways of how 
the three fermions can hop around the hexagon. 

We show now that all twelve ways, which are resulting from 
clockwise/counter-clockwise processes as well as different se- 
quences, have a positive matrix element, (i) Fermion configu- 
rations which result from clockwise or counter-clockwise pro- 
cesses are related by an even number of fermion exchanges, 
e.g. (1,2,3) (2,3,1), and thus have the same sign. In 
passing we mention that this is generally true for ring-hopping 
processes which involve an odd number of fermions, while 
processes with an even number of fermions have the opposite 
sign and thus cancel each other, (ii) the matrix-element does 
not depend on the sequence in which the fermions hop. This 
can be seen if we choose "ordered" fermion configurations 
with respect to the tetrahedra, i.e., we fill the vacuum |0) by 

|c)=4,.«---4,a.cLaJO). (A2) 

Because of the hardcore-dimer constraint, there is at most one 
fermion on each tetrahedron i which then has an arbitrary a. 
In Hamiltonian ( lAll i. the fermionic operators come in pairs for 
each i and thus there is always an even number of fermionic 
operators which need to be exchanged to evaluate the matrix 
elements. Thus, the sign does not depend on the sequence in 



which the operator pairs are applied and all non-zero matrix 
elements are equal to +1. The argument works only for the 
dimer model. If we consider a related loop model on the same 
lattice, there are two fermions per tetrahedron (i.e., two dimers 
attached to each site) and the sign problem can no longer be 
avoided using this method. 

Appendix B: Landau theory for R-state 

In this paper we studied the zero-temperature quan- 
tum phase transition from the ordered R-state into a 
quantum [/(l)-liquid. However we note that the finite- 
temperature, phase transitions between ordered and classical 
U{1) "Coulomb" phases in three-dimensional classical dimer 
models have also proved be a very interesting^ii^^. These 
phase transitions do not fall into the conventional Landau 
paradigm, since the high temperature Coulomb phase exhibits 
the power-law decay of correlation functions characteristic of 
a C/(l) gauge theory, rather than the exponential decay pre- 
dicted by a Ginzburg-Landau theory of the low-temperature 
ordered state. 

In this context it is interesting to note that the present order 
parameter symmetry permits a third-order invariant in a Lan - 
dau action, as well as the usual quadratic and quartic terms : 

6 

Ar = a^m|, 

i=l 

+6 (7TiR,iTOR_2mR^3 + mR^4mR,5mR_6) 

6 / 6 \ 2 

+Ci ^ , + C2 X] ) 

i=l \i=l / 

+C3 [(m| i +rnl 2 + m\ r^Y 

, / 2 2,2 2,2 2 \2 
+C4 (mp, ^mp, 4 + JTlR 2™R,5 + "^R.gmp g) 

(Bl) 

In the ordered phase with (tor) = m ^ 0, the cubic term 
selects eight possible configurations. For 6 < these have 
the form given in Table HV] and are of course the eight R-state 
configurations, which become exact ground states of the quan- 
tum problem in the limit /i ^ — oo. For 6 > the cubic term 
does not select a valid (static) dimer configuration. 

Appendix C: Comparison of variational wave functions 

In this Appendix we explore the transition from solid to liq- 
uid phases of the quantum dimer model on a diamond lattice 
using a variety of variational wave functions. This is an es- 
sential preliminary to Green's function Monte Carlo (GFMC) 
simulations, which rely on importance sampling of configura- 
tions to achieve convergence in a finite time. Except where 
stated otherwise, the calculations presented here were per- 
formed using variational Monte Carlo, within the stochastic 
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TABLE IV: Eight R-states selected by the cubic invariant 6 < in 
Eq. l IBlb . expressed in terms of the coefficients of the order parameter 
Eq. l llSt , within an ordered state with m — (mR/\/3)- The four 
R-states labelled = 1, . . . , 4 have the opposite chirality from those 
labelled C = 5, . . . , 8. 
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reconfiguration algorith m^ for a 1024-bond cluster with 
cubic symmetry. 




FIG. 20: (Color online) Value of jj, for which the one-parameter 
variational wave function [cf. Eq. JCll ll has the lowest energy, 

as a function of the variational parameter a. Results are plotted for 
cubic clusters ranging in size from 128 to 2000 bonds. 



Perhaps the simplest wave function we can consider is 

\ipa^')<p,c = exp[aJ\ff] |V'Rk)0.c 



where a is the single variational parameter associated with 
the total number of flippable hexagons A/f, and \ipRK)<t,,c is 
the RK point ground state function, i.e. the equally weighted 
sum of all dimer configurations within a given subsector of 
the Hilbert space. In the discussion of ground state proper- 
ties which follows, we will consider the maximally connected 
subsector for </> = 0, which comprises all dimer configurations 
connected to an R-state by matrix elements of Eq. ([T]). 

We can view Eq. dCll ). as the Jastrow form which follows 
from Eq. (l22T l. For this reason we anticipate that this wave 
function will work well in the vicinity of the RK point, but will 
tend to over-emphasize any liquid phase occurring for /i 1. 
The wave function iV'^^O^.c h^s the advantage that the value 



of f.1 which minimizes the ground state energy for a given a 
can be found by evaluating averages at the RK point 



i(M/Cf + /CfA/'f)RK - (/Cf)RK(A/'f) 



RK 



(Af )rk - yvt/RK 



mi 



(C2) 



where /Cf is defined by Eq. (|2|. In Fig.|20]we present results 
for /i as a function of a for cubic clusters ranging in size from 
128 to 1024 diamond lattice bonds. It is tempting to interpret 
the abrupt change in dfi/da at a « 0.4, ^ w —0.5 as evi- 
dence of a first order phase transition. However comparison 
with exact diagonalization, and with other variational wave 
functions below, suggests that gives very poor energies 

away from the RK point, and is only really useful for /i ^ 1. 
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FIG. 21: (Color online) Comparison of different variational wave 
functions for 1024-bond cluster: (a) order parameter and (b) energy 
per site: (green) It/'q"') [cf. Eq. (|UD] (orange) [cf. Eq. Oil: 

^Cl) (red) i^^^;)[cf. Eq. d); (black) jV'X) [^f- Eq. Ml 



We can construct another simple one-parameter variational 
wave function based on the order parameter jtir for the R-state 



(C3) 



This wave function gives a good account of the ordered R- 
state, and exhibits a clear first order phase transition transition 
from as a jump in (ttir) for ^ sa 0.25 in a 1024-bond cubic 
cluster — cf. Fig.|2lla). However from Fig.llHb) we can see 
that IV'^^'^) has a lower energy approaching the RK point. 

Seeking the advantages of both terms, we try a two- 
parameter wave function 



(C4) 
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FIG. 22: (Color online) Quantum hysteresis — variational Monte 
Carlo results for IV'a'^'^) [cf. Eq. ^] in a 1024-bond cluster, (left) 
energy per site as a function of the order parameter rriR, (right) order 
parameter as a function of /i. Simulations converge to different local 
minima for a range of values of 0.28 < /i < 0.43, according to 
whether they are started from an ordered or a disordered state. 



FIG. 23: (Color online) Energy per site as a function of fi. GFMC 
results for 1024-bond cluster (points) are compared to Fade approx- 
imants for series expansion around an ordered state (solid lines) and 
perturbation expansion in the vicinity of the RK point (dashed lines). 



This gives lower energies overall, and predicts a transition 
from solid to liquid for fi ss 0.2 in a (1024 bonds). However 
this wave function does a poor job of describing the quantum 
melting of the R-state. 

In order to better model this, we introduce correlations be- 
tween individual flippable plaquettes, on which all dynamics 
in the QDM depend 



exp[aA/f + f^niR 



E 



(C5) 



Here = 1 if the plaquette i is flippable and = 0, oth- 
erwise. For complete generality we should also allow for the 
chirality of the flippable plaquette, and consider a separate 7^ 
for each indistinguishable pair of plaquettes in the cluster. For 
the smallest 128 bond cluster there are 16 of these, and the 
number grows to 84 in the 1024-bond cluster However in 
practice we do not find that including chirality brings signifi- 
cant gains in energy, and we restrict our calculations to a set 
of about 40 7y . 

This wave function has a significantly lower energy than 
any of the alternatives considered above for —0.5 < < 0.25, 
and predicts a transition for ^ « 0.4 (1024-bond cluster — 
black lines in Fig. I2TI ). The wave function also treats fluctu- 
ations accurately enough to capture some a part of the hys- 
teresis associated the quantum melting of a solid. In Fig. |22] 
we show the order parameter and energy for the 1024-bond 
cluster, obtained using iV'a^'^)' seeding VMC from different 
states. Starting from an ordered state, we obtain an ordered 
state solution for ^ < 0.42, while starting from a disordered 
configuration we could see a minimum for /i > 0.3. These re- 
sults confirm the first order character of the phase transition, 
and extend the usefulness of IV'q^'^) a guide function for 
GFMC some way into the region where it predicts a liquid 
phase. 




FIG. 24: (Color online) Monopole string tension ki = Ai/L, plot- 
ted for a range of /i spanning the ordered R-state phase. Solid lines 
denote different Fade approximants to the series expansion about an 
R-state configuration [cf. Eq. ( ID2H . Foints denote the results of 
GFMC simulations of 432-bond (red) and 1024-bond (green) cubic 
clusters. 



Appendix D: Series expansion about the ordered R-state 



There are two limits in which the exact ground state of 
Eq. ([T]) is known. One is the RK point fi — 1, where the 
ground state is the equally weighted sum over all dimer con- 
figurations. The other is ^ — !■ —00, where the ground state is 
one of eight degenerate R-state configurations. In the second 
case we can construct a perturbation theory in by expand- 
ing fluctuations about an ordered R-state. We find that the 
ground state energy (per bond) in the zero-flux sector behaves 
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The first term in this expansion counts the number of flippable 
plaquettes in the R-state (a quarter of the number of bonds). 
Subsequent terms are derived combinatorially, by counting the 
number of possible ways of returning to the the initial R-state, 
using a given number of applications of /Cf and following the 
recipe of the Rayleigh-Schodinger perturbation theory. This 
can be calculated for a finite-size cluster, or an infinite lattice; 
the results quoted here are for an infinite lattice. The resulting 
series Eq. (ID 11 1 contains only terms of odd order in l//i. 

The series for /N [Eq. (IDlb 1 can be continued to values 
of n relevant to our simulations by means of Pade approxi- 
mants. In Fig.|23]we compare several different Pade approx- 
imants to Eq. (IDlb . with the ground state energy per bond 
found in GFMC simulations of a 1024-bond cluster We find 
essentially perfect agreement between our simulations results 



and the results of perturbation theory for /i < —0.8, at which 
point the different Pade approximants start to diverge. Some 
of the Pade approximants remain in agreement with the sim- 
ulations up to i^L w —0.5. For comparison, we also include 
results for Eq/N derived in an expansion about the RX point, 
as described in Section llll Dl Again, the agreement between 
perturbation theory and simulation is found to be very good in 
the limit /i — > 1. 

For /i — > — oo, it is also possible to perform a series about 
a state with unit-flux, i.e. an R-state with a single string exci- 
tation. From Eq. (l2Tl i we find the monopole string tension to 
be 



Kl 



-2/i 



1 1 
15/1 



71243 1 



1995840 /i3 
1044561213049 1 
"41535007113600000/1^ ^ 



(D2) 



In Fig. |24] we compare the perturbational expansion and 
GFMC simulation results for the string tension ki. The pre- 
diction of the perturbation theory for /i < — 1 matches well to 
the results of simulations for /i > — 1. 
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We note that these, arbitrary, units differ from those defined in 
Eq. ([HI* by a factor 4/\/3. 



